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Study of digital pulse shape discrimination method for n -7 separation 

of EJ-301 liquid scintillation detector* * 
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Abstract: A digital pulse shape discrimination system based on a programmable module NI-5772 has been estab¬ 
lished and tested with EJ-301 liquid scintillation detector. The module was operated by means of running programs 
developed in LabVIEW with the sampling frequency up to 1.6GS/s. Standard gamma sources ^^Na, ^®^Cs and ®°Co 
were used to calibrate the EJ-301 liquid scintillation detector, and the gamma response function has been obtained. 
Digital algorithms for charge comparison method and zero-crossing method have been developed. The experimental 
results showed that both digital signal processing (DSP) algorithms could discriminate neutrons from gamma-rays. 
Moreover, the zero-crossing method shows better n-7 discrimination at 80 keVee and lower, whereas the charge com¬ 
parison method gives better results at higher thresholds. In addition, the figure-of-merit (FOM) of two different 
dimension detectors were extracted at 9 energy thresholds, and it was found that the smaller one presented a better 
n-7 separation property for fission neutrons. 
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1 Introduction 

For the last decades, concepts of accelerator driving 
subcritical systems (ADS) and spallation neutron sources 
(SNS) have been pressed and constructed based on the 
spallation reactions [ll-l3|. High intensity neutron flux 
and gamma-rays would be produced around the spalla¬ 
tion targets. Kinds of neutron detectors have been used 
to monitor the leakage neutrons in real time in mixed 
neutrons-gamma-rays radiation environments, such as 
fission chamber, ®He counter and organic liquid scintilla¬ 
tion detector et al. Among of these detectors, the EJ-301 
liquid scintillation detectors have been widely employed 
because of their excellent neutron-gamma discrimina¬ 
tion, high efficiency for fast neutron detection and the 
superior time resolution. Aiming to discriminate neu¬ 
trons from gamma-rays for the EJ-301 liquid scintillation 
detector, two most common pulse shape discrimination 
(PSD) methods could be used. One of them is the zero¬ 
crossing method, which extracts the zero-crossing time 
of the doubly differential bipolar pulse Another 


is the charge comparison method based on independent 
measurements of the integrated charge over two different 
time regions of the pulse^, 0 - 

Recently, with the development of the field pro¬ 
grammable gate array (FPGA) technology and computer 
GPUs, the digital signal processing (DSP) is now possi¬ 
ble. The major difference between analog and DSP tech¬ 
niques is that with the digital method, the current pulse 
from the anode of photomultiplier tube (PMT) is digi¬ 
tized immediately and all operations are carried out in 
software package. For the analog PSD methods, the op¬ 
erations have to be implemented based on CAMAC or 
VME modules along with a series of complicated ana¬ 
log circuits. The DSP system offers significant advan¬ 
tages over the analog system in the area of conveniences, 
real-time properties and cramped construction by elimi¬ 
nating extra electronic modules. Moreover, prior to this 
work, some research groups have investigated the DSP 
method, and performed the comparisons of the n -7 sepa¬ 
ration results using the zero-crossiM method and charge 
comparison method, respectively [al^. However, their 
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experimental results showed different behaviours on the 
ability of n -7 discrimination. For example, when com¬ 
paring the zero-crossing method and charge comparison 
method, results in Ref. @ suggested the former showed 
better n^ separation property, whereas the conclusion 
in Ref. [9| was contrary. Therefore, it is very essential 
to develop DSP method, and timely to reinvestigate the 
zero-crossing method and charge comparison method. 

In this paper, a digital acquisition system based on 
the NI-5772 adapter module has been developed. The al¬ 
gorithms of digital pulse shape discrimination method for 
n -7 separation of EJ-301 liquid scintillation detector were 
described in Section 2. The liquid scintillation detector 
was calibrated by the standard gamma-rays sources, and 
the capacity of n -7 discrimination for the new system 
was tested with 252 qj neutron source. At the same time, 
the capacity of n -7 discrimination for two different di¬ 
mension liquid scintillation detectors was also studied. 


2 Experimental details 


2.1 Algorithm for PSD 


Fig.UDa) presents the typical current pulses related to 
neutrons and gamma-rays. Pulse generated by gamma- 
ray decays faster to the baseline than the neutron in¬ 
duced pulse. The major difference between these two 
pulses occurs in their tail. It is this fact that neutrons 
and gamma-rays could be discriminated by analyzing the 
contributions of slow components to the total light out¬ 
put. 

In this paper, two types of digital data processing 
techniques are employed for n -7 discrimination. One of 
them is the digital implementation of the conventional 
charge comparison method. Qtotai and Qsiow represent 
the total charge of the current pulse and the charge in 
the slow components, respectively, as shown in Fig.[T])a). 
Thus neutrons and gamma-rays could be separated ac¬ 
curately in mixed radiation fields through comparing the 
difference of Qsiow for each pulse. The other algorithm is 
zero-crossing method. In Fig. [lib), the bipolar pulse is 
obtained by analyzing a digital differentiator-integrator- 
integrator (CiRi-(R 2 C 2 )^) algorithm to the PMT cur¬ 
rent signal. Therefore, the difference in the tail of dif¬ 
ferent PMT signals is reflected in the zero-crossing time 
of the bipolar pulses. To exploit the difference in the 
zero-crossing time as a parameter for n -7 discrimina¬ 
tion, a digital algorithm acting as a constant fraction dis¬ 
criminator (CFD) is used to determine the zero-crossing 
points. 



Time (0.6ns/channel) 


7 



Fig. 1. (a) Typical PMT signals induced by neu¬ 

trons and gamma-rays. The gamma-ray pulse de¬ 
cays more quickly than the neutron pulse, thus 
small difference can be seen in the tail, (b) PMT 
signals after the pulse shaping process. The 7 -ray 
and neutron pulses cross the zero line at different 
times. 

2.2 Experiment setup 



PXI 

Fig. 2. The diagram of experiment setup. 


The diagram of experimental setup and electronics 
are shown in Fig. |21 A cylindrical EJ-301 liquid scintil¬ 
lation detector of 2 -inch in diameter and 2 -inch long cou¬ 
pled to a ET 9813KB PMT was used to perform the mea¬ 
surements. The detector was operated with a negative 
voltage of 1700V. The current pulses from anode of PMT 
were transferred to an attenuator and then directly dig¬ 
itized using a digital oscilloscope (Nl-5772), which could 
digitize the waveform by running dedicated Lab VIEW 
data acquisition package. Since this acquisition system 
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could only digitize the signals whose maximum voltage 
range was between -IV and IV, the 9dB attenuator was 
employed. In the present work, the digital oscilloscope 
worked at self-triggering mode with a sampling rate of 
1.6GS/s and 12 bit resolution. A neutron source 
with the intensity of 1 . 0 x 10 "* n/s was used to test the n- 
7 discrimination capability of this system, and standard 
gamma sources were also used for energy calibration of 
EJ-301 liquid scintillation detector. Off-line data analy¬ 
sis was performed by running ROOT script files written 
in C-I--I- language. In addition, to compare the capability 
of n -7 discrimination for scintillator in different volume, 
another larger liquid scintillation detector of 2 -inch in 
diameter and 4-inch long was also employed. 

2.3 Energy calibration 


In order to know the specific relationship between 
light output distribution and the realistic energy de¬ 
posited in scintillator, and determine the effective neu¬ 
tron detection thresholds, the gamma energy calibration 
for the EJ-301 liquid scintillation detector was done with 
standard gamma sources ^^Na, *®*’Cs and ®“Co. Because 
the elements of liquid scintillator have very low photo¬ 
electric effect cross section, the obtained spectrum are 
mainly due to Compton scattering electrons 0 - The 
position of the maximum of Compton scattering elec¬ 
trons could be accurately determined through comparing 
the measured light output spectrum with simulated ones 
using different Monte Carlo simulation packages, such as 
CEANT4 [iil,[i3, FLUKA O, Q and CRESP7 [H]. 
Among of these packages, it has been proved that the 
CRESP7 code is very reliable for ener gy calibration. 
More details have been described in Ref. 13 US- 

In this work, the light output spectrum of the detec¬ 
tor for ^^Na, *®*’Cs and ®'*Co gamma sources were simu¬ 
lated by GRESP7 code. Fig. [3] shows the experimental 
and GRESP7 simulated Compton electron spectrum of 
*®*'Cs source. As shown in Fig.|31 the maximum Compton 
edge is corresponding to 81% of the Compton electron 
distribution. The results of calibration together with 
the light output of ^^Na, ^^’^Cs and ®°Co are shown in 
Fig.ia in which the experimental data is fitted by linear 
polynomial. 

The gamma response function of EJ-301 liquid scin¬ 
tillation detector is expressed by: 


L = 4.94 0.095, (1) 

where L is the light output, is the deposited electron 
energy for liquid scintillator in MeV. This calibration is 
implemented to calculate the equivalent electron energy 
(MeVee), where 1 MeVee corresponds to the total light 
output induced by a 1 MeV electron. 



Channel number 


Fig. 3. Simulated Compton electron spectrum of 
*^*"08 source (blue dash line) and the same spec¬ 
trum folded with energy resolution (red solid 
line). Green solid line is experimental light output 
of «^Cs. 



Qtotal 


Fig. 4. Experimental light output for ^^Na, *®*^Cs 
and ®°Co gamma sources. The inset shows the 
calibration result using the Compton edges from 
the detected spectrum. 


3 Results and discussion 

3.1 Charge comparison method 

Two-dimensional scatter plot of the n -7 discrimi¬ 
nation using charge comparison method at an energy 
threshold of 180 keVee is shown in Fig. [5l There are 
two separated strips in the 2D graph where the upper 
one is connected with neutrons because of larger Qsio-w 
for equal light output. The overlap is due to the case 
of very low energy neutrons and gamma rays, which is 
the region that determines how well the neutrons are 
discriminated from gamma-rays. 
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Fig. 5. 2D plot of Qtotai versus Qsiow for 
source, the energy threshold is 180 keVee. 


The quality of n -7 discrimination depend on several 
reasons, such as the employed algorithms, the sampling 
rate and the type of detectors and so on. In order to 
check the character of n -7 separation, a figure-of-merit 
(FOM) is introduced: 


FOM = 


A 

A 7 + An ’ 


( 2 ) 


where A is the separation between the peaks of the neu¬ 
tron and gamma events. The A 7 and An are the FWHM 
(full width at half maximum) of the gamma and neu¬ 
tron peaks, respectively. The procedure of extracting the 
FOM is illustrated with ID histogram of Qsiow/Qtotai, as 
shown in Fig. | 6 l 



Qsiow/Qtotai 


Fig. 6. The procedure of calculating the FOM. 

FOMs over different time intervals were extracted for 
the purpose of evaluating the optimum integration in¬ 
tervals for Qsiow and Qtotai- The integration intervals 
were optimized for minimizing the number of misclassi- 
fied pulses. In this work, the total integration was cal¬ 
culated from beginning of the pulse to the 50ns after the 
pulse maximum for all the signals. Meanwhile, FOMs 
for different time intervals of Qsiow are shown in Fig. [7] 


at 5 energy thresholds (dOkeVee, 120keVee, ISOkeVee, 
240keVee and SOOkeVee). It is found that the optimum 
n -7 separation could be derived since the tail integration 
started 21ns and ended 50ns from the pulse maximum 
for the all studied energy thresholds. The extracted op¬ 
timum integration intervals will be used for charge com¬ 
parison algorithm in following parts. 



Fig. 7. FOMs for different integration intervals 
of Qsiow, the energy thresholds are GOkeVee, 
120keVee, ISOkeVee, 240keVee and SOOkeVee. 

3.2 Zero-crossing method 



Qtota! 


Fig. 8. 2D plot of zero-crossing time versus Qtotai 
for ^®^Cf source, the energy threshold is 180 
keVee. 


A series of FOMs were calculated at different con¬ 
stants for CiRi and R 2 C 2 , and it was found that the 
optimum n -7 discrimination was achieved when the dif¬ 
ferentiation time constant (CiRi) was 5ns and the inte¬ 
gration time constant (R 2 C 2 ) was 30ns. Fig. |5] exhibits 
the zero-crossing time distribution versus Qtotai for the 
same experimental data set. The upper peak corresponds 
to neutrons because the neutron pulse crosses the base¬ 
line much later than the gamma pulse. The FOM in 
this case is 0.91. In order to investigate the effect of en¬ 
ergy threshold on the quality of n -7 separation, FOMs at 
11 different energy thresholds have been calculated and 
shown in Fig. |9l The FOMs are improved from 0.78 to 
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1.04 when the thresholds increase from 30 keVee to 600 
keVee. This phenomenon could be explained that the 
influence of electronic noise decreases with larger energy 
threshold. 

FOMs calculated using charge comparison algorithm 
for the same energy thresholds were also shown in Fig. |9] 
so as to compare the n-y separation capacity of charge 
comparison and zero-crossing technique. The contrastive 
results show that the zero-crossing method is better for 
energy thresholds of 80 keVee and lower. As for thresh¬ 
olds higher than SOkeVee, the charge comparison method 
exhibits optimal n-y separation property. Therefore, the 
zero-crossing method is suitable for outputs lower than 
80 keVee and charge comparison method is the best 
choice for higher equivalent electron energy. The re¬ 
sults are similar with Ref. except that they inves¬ 
tigated the comparison between charge integration and 
zero-crossing method using analog PSD techniques. 



Fig. 9. Comparison of FOMs between charge inte¬ 
gration and zero-crossing method at differen en¬ 
ergy thresholds. 


3.3 PSD for different sized detectors 

In addition, two liquid scintillation detectors men¬ 
tioned in section 2.2 were tested using source to 

compare the capacity of n-y discrimination for detectors 
in different size. The contrastive results calculated us¬ 
ing charge comparison algorithm are shown in Fig. 1101 
Significantly, the smaller liquid scintillation detector has 
greater FOM than the larger one at every energy thresh¬ 
old. This is perhaps because of longer path of photon 
transportation in liquid scintillation for the larger one, 
which distorts the time information at the tail of pulses. 
Consequently, for fission neutrons, the small dimension 
detectors are more universally used to minimize the dis¬ 
tortion of neutron properties. 



Threshold (keVee) 

Fig. 10. Comparison of FOM for different size de¬ 
tectors at various energy thresholds. 

4 Summary 

A digital acquisition system based on the NI-5772 
adapter module for n-y discrimination was established 
and tested with a ^®^Cf neutron source. This system 
eliminates the need for QDC, TAG and delay cable ect. 
used in analog PSD techniques. The integration charge 
and zero-crossing time could be extracted through pro¬ 
cessing the digital waveforms off-line. 

The energy calibration for EJ-301 liquid scintillation 
detector was done via ^^Na, ^^’^Cs and ®“Co sources, 
and the gamma response function was obtained as: 
L=4.94*i?e-0.095. Two different digital PSD algorithms 
were used to perform the n-y discrimination of EJ- 
301 liquid scintillation detector: (a) charge compari¬ 
son method (b) zero-crossing method. Both algorithms 
clearly showed the power of a digital system in achiev¬ 
ing good PSD. The digital charge comparison method 
presented the optimum n-y separation when the inte¬ 
gration of slow components started 21ns after the pulse 
maximum. At energy thresholds 80 keVee and lower, 
the zero-crossing method gave better n-y discrimination. 
At higher energy, the charge comparison method pre¬ 
sented better separation between neutron and gamma- 
ray events. In addition, as for liquid scintillation detec¬ 
tors in different dimension, the smaller one showed bet¬ 
ter n-y discrimination property for fission neutrons. The 
experimental results showed that such digital signal pro¬ 
cessing technique could be used for neutron monitoring 
of ADS and SNS in the future. 
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